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Affine equivariant rank-weighted L-estimation 
of multivariate location 

Pranab Kumar Sen, Jana Jureckova and Jan Picek 


Abstract In the multivariate one-sample location model, we propose a class of fle¬ 
xible robust, affine-equivariant L-estimators of location, for distributions invoking 
affine-invariance of Mahalanobis distances of individual observations. An involved 
iteration process for their computation is numerically illustrated. 


1 Introduction 


The affine-equivariance and its dual affine-invariance are natural generalizations of 
univariate translation-scale equivariance and invariance notions (Eaton 1983). Con¬ 
sider the group of transformations of Rp to Rp: 

Xh-> Y^BX-fb, bGKp (1) 

where B is a positive definite p x p matrix. Generally with the choice of B we do 
not transform dependent coordinates of X to stochastically independent coordinates 
of Y. This is possible when X has a multi-normal distribution with mean vector 0 
and positive definite dispersion matrix Z, when letting = BB^, so that EY = 
^ = B0 -fb and dispersion matrix BZB^ = Ip- To construct and study the affine 
equivariant estimator of the location 0 we need to consider some affine-invariant 
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(AI) norm. The most well-known affine invariant norm is the Mahalanobis norm, 
whose squared version is 

= {x-d)'^z-\x-e) = \\x-e\\l ( 2 ) 

where E is the dispersion matrix of X. To incorporate this norm, we need to use its 
empirical version based on independent sample Xi,... ,X„, namely 

Sij = ^-{Xt - X,)Tvr ‘ (X,- - X,), !</<;<« 

where V* = {n{n — 1))^* Li<;<j<n(X,' — Xj){Xi — Xj)~^. To avoid redundancy, we 
may consider the reduced set 

4 = (X,-X„)Tv-‘(X,-X„), 1 (3) 

n— 1 

V„= ^(X,--X„)(X,--X„)T (4) 

1=1 

which forms the maximal invariant (MI) with respect to affine transformations ([T]i. 
An equivalent form of the maximal invariant is 

dni={Xi-X„)'^\-\Xi-X„), /= l,...,n (5) 

V„ = f;(X,--X„)(X,--X„)T (6) 

1=1 

(Obenchain (1971)). Note that all the dm are between 0 and 1 and their sum equals to 
p. Because dm are exchangeable, bounded random variables, all nonnegative, with 
a constant sum equal to p, the asymptotic properties of the array {d„i,... ,d„n)~^ 
follow fromChernoff and Teicher (1958) and Weber (1980). Similarly, dm = p- 
Neither X„ nor V„ is robust against outliers and gross errors contamination. As such, 
we are motivated to replace A„ and V„ by suitable robust versions and incorporate 
them in the formulation of a robust affine equivariant estimator of 0. If 0„ is some 
affine equivariant estimator of 0, then writing 

V„ = f;(X,-0„)(X,-0„)T 
1=1 

we may note that V„ is smaller than V„ in the matrix sense. However, it cannot be 
claimed that the Mahalanobis distances (|5]l can be made shorter by using 0„ instead 
of X„, because L/Li (X; — 0„)^ (X, ~ d„) = p. Our motivation is to employ the 

robust Mahalanobis distances in the formulation of robust affine equivariant estima¬ 
tor of 0, through a tricky ranking of the Mahalanobis distances in Q and an iterative 
procedure in updating an affine equivariant robust estimator of 0. 

The robust estimators in the multivariate case, discussed in detail in Jureckova 
et al. (2013), are not automatically affine equivariant. With due emphasize on the 
spatial median, some other estimators were considered by a host of researchers, and 
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we refer to Oja (2010) and Serfling (2010) for a detailed account. Their emphasis is 
on the spatial median and spatial quantile functions defined as follows: 

Let i^p_i(0) be the open unit ball. Then the u-th spatial quantile Qp{\i), u £ 
■^p-i(O) is defined as the solution x = (u) of the equation 



Particularly, (0) is the spatial median. It is equivariant with respect toy^Bx-f 
b, b £ Kp, B positive definite and orthogonal. However, the spatial quantile function 
may not be affine-equivariant for all u. 

Among various approaches to multivariate quantiles we refer to Chakraborty 
(2001), Roelant and van Aelst (2007), Hallin et al. (2010), Kong and Mizera (2012), 
Jureckova et al. (2013), among others. Lopuhaa & Rousseeuw (1991) and Zuo 
(2003, 2004 2006), Lopuhaa & Rousseeuw (1991) and Zuo (2003, 2004 2006), 
among others, studied robust affine-equivariant estimators with high breakdown 
point, based on projection debths. An alternative approach based on the notion of 
the depth function and associated U-statistics has been initiated by Liu at al. (1999). 
Notice that every affine-invariant function of (Xi,... ,X„) depends on the X,- only 
through the maximal invariant; particularly, this applies to the ranks of the dni and 
also to all affine invariant depths considered in the literature. In our formulation, 
affine equivariance property is highlighted and accomplished through a ranking of 
the Mahalanobis distances at various steps. 


2 Affine equivariant linear estimators 

Let X £ Rp be a random vector with a distribution function F. Unless stated other¬ 
wise, we assume that F is absolutely continuous. Consider the qroup of affine 
transformations X i~-> Y = BX + b with B nonsingular of order p x p, b £ Mp. Each 
transformation generates a distribution function G also defined on Rp, which we 
denote G = Fp j,. A vector-valued functional 9{F), designated as a suitable measure 
of location of E, is said to be an affine-equivariant location functional, provided 


9{FB.b) = B0(E) +b Vb £ Rp, B positive definite. 


Let r(E) be a matrix valued functional of F, designated as a measure of the scatter 
of F around its location 0 and capturing its shape in terms of variation and co¬ 
variation of the coordinate variables. r{F) is often termed a covariance functional, 
and a natural requirement is that it is independent of 0(E). It is termed an affine- 
equivariant covariance functional, provided 


r{FB b) = BE(E)B''' Vb £ Rp, B positive definite. 
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We shall construct a class of affine equivariant L-estimators of location pa¬ 
rameter, starting with initial affine-equivariant location estimator and scale func¬ 
tional, and then iterating them to a higher robustness. For simplicity we start 
with the sample mean vector X„ = \ with the matrix y„ = = 

— X„)(X,' — X„)^ = nLn, n > p, where £„ is the sample covariance ma¬ 
trix. Let R„i = Y!j=\t[dni < dni] be the rank of d„i among r/„i,.. i = 1,... ,n, 
and denote R„ = {Rni,- ■ ■ vector of ranks. Because F is continuous, the 

probability of ties is 0, hence the ranks are well defined. Note that dni and R„i 
are affine-invariant, i = Moreover, the R„i are invariant under any strictly 

monotone transformation of dni, / = 1,..., n. Furthermore, each X, is trivially affine- 
equivariant. We introduce the following (Mahalanobis) ordering of Xi,...,X„ : 


X/ ^ Xj dfii ^ dnj, i ^ j — 1,,11. 


(7) 


This affine invariant ordering leads to vector of order statistics X„:i ^ ... ^ X„:„ of 
the sample X„. In the univariate case with the order statistics X„:i < ... < X,,-.,,, we 
can consider the k-order rank weighted mean (Sen (1964)) defined as 


Tnk = 


n 

2k+\ 


-1 n—k 

L 

i=k+l 


I — 

k 



k = 0,.. 


n + 1 
2 


For k = 0, Tnk = Xn and k=\(nF\) /2] leads to the medianX„. In the multivariate 
case, the ordering is induced by the non-negative dnu and the smallest dni corre¬ 
sponds to the smallest outlyingness from the center, or to the nearest neighborhood 
of the center. Keeping that in mind, we can conceive by a sequence {k„} of non¬ 
negative integers, such that is ^ in n, but is \ in n, and for fixed k 

put 

^nk = ( ) E / [Rni < kn] ( ) X,-. 

L„i; is affine-equivariant, because the dni are affine invariant and the X,- are trivially 
affine equivariant. Of particular interest is the case of k = 1, i.e., 

n 

^n\=Ki^Y^i < kn] X/ 

/=1 

representing a trimmed, rank-weighted, nearest neighbor (NN) affine-equivariant 
estimator of 0. In the case k = 2 we have 

Ln2 = ( 2 " ) £ ^ 

which can be rewritten as L „2 = LJLi >VnJ?,„X, with the weight-function 
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Wni = 



{k„-i) ... i= 

... i> k„. 


We see that L„ puts greater influence for /?,„■ = 1 or 2, and Wnk„ = 0; w„i = 2/A:„. 
For k>3, even greater weights will be given to /?„, = 1 or 2, etc. For large n we can 
use the Poisson weights, following Chaubey and Sen (1996): 

A typical X is chosen somewhere in the middle of [0,1]. Then L[] = Yfi=i 
represents an untrimmed smooth affine-equivariant L-estimator of 0; for A 0 we 
get the median affine-equivariant estimator, while A 1 gives a version of L„ 2 -type 
estimator. If A is chosen close to 1/2 and = o{y/n), then tail L;>yt„ converges 
exponentially to 0, implying a fast negligibility of the tail. Parallelly, the weights 
w„{i) can be chosen as the nonincreasing rank scores a„(l) > a„(2) > ... > a„{n). 

To diminish the influence of the initial estimators, we can recursively continue in 
the same way: Put = L„ and define in the next step: 


1=1 

41’ = (X,--L,4)T(A,4)^'(X,--Li'’ 

7(1)1 


(i)l 


(O'lT 


;=1 


The second-step estimator is = Yli=\ proceed, so at 

the r-th step we define A„ , , 1 <i <n and the ranks R„ ' analogously, and get 

the r-step estimator 

Li^’=Ew„(4r‘’)Xo ^>1- (8) 

i=l 

(r) 

Note that the are affine-invariant for every 1 < / < n and for every r > 0. Hence, 
applying an affine transformation Y, = BX, -fb, b £ Rp, B positive definite, we see 
that 

lW (Yi ,..., Y„) = BlW (Xi,...,X„) + b. (9) 

Hence, the estimating procedure preserves the affine equivariance at each step and 

(r) 

L„ ' is an affine-equivariant L-estimator of 0 for every r. 


The algorithm proceeds as follows: 

(1) Calculate X„ and Ai°’ = 1 (X,- - X„) (X; - X„) ^. 

(2) Calculate 4-’= (X,-X„)T(A<i'’’)-i(Xi-X„), !</<«. 
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(3) Determine the rank of among ,...,dm , i = 1, ■ ■ ■ ,n. 

(4) Calculate the scores a„{i), i = 1,... ,n 

(5) Calculate the first-step estimator 

(6) Ai') 

(7) 4'^ = (X,'-LWr(Al‘>r'(X,--LW), l<i<n. 

(8) = the rank of d^]^ among d^]^di]}, i=l,...,n. 

(9) 

(10) Repeat the steps (6)-(9). 

(r) 

The estimator L„ ' is a linear combination of order statistics corresponding to 
independent random vectors Xi,..., X„, with random coefficients based on the ex- 

(r) (r) 

changeable . The asymptotic distribution of L„ ' under fixed r and for « oo is a 

(r) 

problem for a future study, along with the asymptotic distribution of the d^/ and of 
the rank statistics. For the moment, let us briefly recapitulate some asymptotic prop- 

(r) (r) fr) 

erties of the . Note that = p V r > 0, and that the d^f are exchangeable 

nonnegative random variables with a constant sum and ~ n every r>0. 

Let and S* = (X, - 0)^£-^ (X,-- 0), l<i< 

II, r > 1 where E is the covariance matrix of Xi. Let G„ Ty) = P{nd),^ < y} be 
the distribution function of the nd\!^} and letGi'^^(y) = «^* y S 

be the empirical distribution function. Side by side, let G^iy) = P{^ni — y} 

G* (y) = P{5* < y} be the distribution function of and 5* respectively. By the 
Slutzky theorem, 

|G:,(y)-G*(y)|^0 as « ^ oo. 

Moreover, by the Courant theorem. 


TAx 


Chmin(AB ^) = inf < sup 


Bx 


x^Ax 

x^Bx 


:Ch^ax(AB-'), 


we have 


max 

\<i<n 


ndl!) 


5^ 


-1 


< max 


Chn 


-( A ^‘0 


Wt-i^ 


£]-l 


Ch„ 




-1 


so that 


iAi' 

n 


(r) 


|gW-g: 


0. In a similar way, 


0||, where the right-hand side is Op{rr^/^). Because d\^) are exchangeable, bounded 
and nonnegative random variables, one can use the Hoeffding (1963) inequality to 
verify that for every Cn > 0, there exist positive constants K\ and v for which 


-vnc2 


^{|Gl''^(y)-Gl''^(y)l>c„}<^fie 
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Thus, using c„ = O {n^/^logn) can make the right-hand side to converge at any 
power rate with This leads to the following lemma. 

Lemma 1. As fl —y CXD, 

suPrfeR+ - Gi"\y) - Gi[\y') + Gi’'\y')\ : ly-y'l < n^'/V21ogn| 

“= 0{n^^logn). (10) 

Proof (outline). The lemma follows from the BoreTCantelli lemma, when we notice 
that both (y) and Gn^ (y) are /^inyG R+, and that (0) = gJ,^^ (0), Gi’’'^ (°°) = 
GW(oo) = 1. 

Theorem 1. Let 


Wnr{t) = 1(f)) _ f g [0^ 1], ^ l^^^(f). 0 < f < 1}. 


Then Wnr W in the Skorokhod ^[0, 1] topology, where W is a Brownian Bridge 
on [0,1]. 

Proof (outline). The tightness part of the proof follows from Lemma[T] For the con¬ 
vergence of finite-dimensional distributions, we appeal to the central limit theorem 
for interchangeable random variables of Chernoff and Teicher (1958). 


If the X; have multinormal distribution, then 5* has the Chi-squared distribution 
with p degrees of freedom. If X, is elliptically symmetric, then its density depends 
on h{\\x - ellx), with h{y), y>Q depending only on the norm ||y||. If p > 2, as it is 
in our case, it may be reasonable to assume that H{y) = Jq h(u)du behaves as ~ yP!'^ 
(or higher power) for y —0. Thus y ~ [//(y )]^/p (or [H{y)Y, r < 2/p) for y —0. 
On the other hand, since our choice is k„ = o{n), the proposed estimators L„i and 
L „2 both depend on the X, with dni of lower rank {R„i < k„ or n^^Rni < n^^k„ 0). 

Hence, both L„i and L „2 are close to the induced vector Xpj where [1] = {/: /?„, = 
1}. If the initial estimator is chosen as Xpj and A„[i] = n^^Y!i=i0^i ~ X[i])(X, — 
Xp])^, then the iteration process will be comparatively faster than if we start with 
the initial estimators X„ and 

The proposed L„i, L „2 are both affine equivariant and robust. If we define the 
D-efficiency 


D 


('-) _ 



!•< 1, 


( 11 ) 


then it will be slightly better than that of the spatial median; the classical X„ has the 
best efficiency for multinormal distribution but it is much less robust than L„i and 

L^2- 
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3 Numerical illustration 


3.1 Multivariate normal distribution 


The procedure is illustrated on samples of size u = 100 simulated from the normal 
distribution ^(0,Z) with with 


0 = 



L = 


1 Xjl Ijl 
1/2 1 1/2 
1 / 21/2 1 


( 12 ) 


and each time the affine-equivariant trimmed L„i-estimator (A:„ = 15) and affine- 
equivariant L„ 2 -estimator were calculated in 10 iterations of the initial estimator. 
5 000 replications of the model were simulated and also the mean was computed, 
for the sake of comparison. Results are summarized in Table 1. Figure 1 illustrates 
the distribution of estimated parameters 0i, 02 , 03 for various iterations of L„i- 
estimator and L„ 2 -estimator and compares them with the mean and median. Tables 
2-3 and Figure 2 compare the D-efficiency of proposed estimators. 

The Mahalanobis distance is also illustrated. One sample of size n = 100 was 
simulated from the bivariate normal distribution with the above parameters. After¬ 
wards, the Mahalanobis distances di\ = (Xj — XY {Xj — X), i = were 

calculated. They represent n co-axial ellipses centered at X - see Figure 3 (black 
ellipses). The modified Mahalanobis distances replaced X by the affine-equivariant 
trimmed L„i-estimator with = 15 (see the blue ellipses on Figure 3) and affine- 
equivariant L„ 2 -estimator (see the red ellipses on Figure 3) with analogous modifi¬ 
cation of S,f 


Table 1. Normal distribution: The mean in the sample of5 000 replications of 
estimators L„i (trimmed) andL,„ 2 , sample sizes n = 100 


i 

^(1 

*^n2 

1 

0.999607 2.001435 -0.998297 

0.999401 1.999796-0.999943 

2 

0.999473 2.001423 -0.996185 

0.999083 1.999584 -0.999782 

3 

0.999519 2.000290 -0.993274 

0.998926 1.999509 -0.999801 

4 

0.999435 2.000190 -0.991901 

0.998854 1.999496 -0.999871 

5 

0.999474 2.000771 -0.991295 

0.998811 1.999476-0.999973 

6 

0.999646 2.001285 -0.990964 

0.998781 1.999471 -1.000032 

7 

0.999926 2.001519 -0.990829 

0.998773 1.999470-1.000049 

8 

0.999952 2.001529 -0.990803 

0.998772 1.999472 -1.000068 

9 

0.999939 2.001497 -0.990745 

0.998775 1.999489-1.000071 

10 

0.999853 2.001424 -0.990711 

0.998779 1.999497 -1.000061 
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Table 2. Normal distribution: The mean, median and minimum of D-efficiency in 
the sample of5 000 replications of estimators L„i (trimmed) andT,,,!, sample sizes 

n = 100 



^n\ 


iterration 

mean median minimum 

mean median minimum 

2 

1.000637 0.999615 0.871029 

0.999862 0.999670 0.951057 

3 

1.001272 0.998593 0.818591 

0.999625 0.999451 0.946809 

4 

1.001106 0.997822 0.793030 

0.999392 0.999231 0.946575 

5 

1.000697 0.997283 0.794157 

0.999205 0.999041 0.946332 

6 

1.000394 0.997211 0.793903 

0.999078 0.998886 0.946022 

7 

1.000131 0.997122 0.793903 

0.998997 0.998802 0.945613 

8 

0.999916 0.996964 0.793903 

0.998951 0.998807 0.945032 

9 

0.999882 0.996924 0.793903 

0.998921 0.998768 0.944518 

10 

0.999860 0.996924 0.793903 

0.998899 0.998706 0.944272 


Table 3. Normal distribution: The 25%-quantile, 75%-quantile and max of 
D-efficiency in the sample of5000 replications of estimators L„i (trimmed) and 
L„ 2 , sample sizes « = 100 



Hi 

H2 

iteration 

25%-quantile 75%-quantile 

max 

25%-quantile 75%-quantile 

max 

2 

0.971669 

1.028444 

1.169137 

0.990914 

1.008819 

1.057115 

3 

0.960891 

1.039437 

1.295245 

0.989675 

1.009489 

1.065753 

4 

0.956068 

1.042774 

1.294952 

0.989358 

1.009426 

1.068542 

5 

0.953186 

1.043542 

1.301147 

0.989098 

1.009239 

1.069236 

6 

0.952193 

1.044435 

1.305327 

0.989000 

1.009182 

1.069437 

7 

0.951535 

1.044942 

1.326996 

0.988916 

1.009142 

1.069600 

8 

0.951441 

1.044394 

1.330791 

0.988839 

1.009109 

1.069226 

9 

0.951452 

1.044562 

1.330791 

0.988802 

1.009087 

1.069236 

10 

0.951356 

1.044749 

1.330791 

0.988771 

1.009062 

1.069278 
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3.2 Multivariate t-distribution 

Similarly, we illustrate the procedure on samples of size n = 100 simulated from 
the multivariate t distribution with 3 degree of freedom with the same 

parameters as in (O. Each time, 10 iterations of affine-equivariant trimmed L„i- 
estimator {kn = 15) and of affine-equivariant L„ 2 -estimator, started from the initial 
estimator, were calculated. 5 000 replications of the model were simulated and the 
mean was computed, for the sake of comparison. Results are summarized in Table 
4. Figure 4 illustrates the distribution of estimated parameters 0i, O 2 , 63 for various 
iterations of L„i-estimator and L„ 2 -estimator and compares them with the mean and 
median. The Tables 5-6 and Figure 5 compare the D-efficiencies of the proposed 
estimators and Figure 6 illustrates the Mahalanobis distances. 


Table 4. t-distribution: The mean in the sample of 5 000 replications of estimators 
L„i (trimmed) and L„ 2 , sample sizes n = 100 


i 

tB 

^n\ 

*^n2 

1 

1.004760 2.002252 -0.992618 

1.003081 2.001199 -0.998781 

2 

1.005034 2.001851 -0.991150 

1.002154 2.000044 -0.999996 

3 

1.005473 2.001430 -0.991330 

1.001744 1.999661 -1.000627 

4 

1.005334 2.000732 -0.992535 

1.001594 1.999526-1.000909 

5 

1.005351 2.000951 -0.993532 

1.001540 1.999475 -1.001028 

6 

1.005007 2.001127 -0.994361 

1.001528 1.999452-1.001069 

7 

1.004636 2.001009 -0.994817 

1.001522 1.999447 -1.001086 

8 

1.004520 2.000930 -0.995011 

1.001524 1.999443 -1.001090 

9 

1.004466 2.000911 -0.995172 

1.001526 1.999442 -1.001090 

10 

1.004445 2.000821 -0.995221 

1.001527 1.999444-1.001086 


Table 5. t-distribution: The mean, median and minimum of D-efficiency in the 
sample of5 000 replications of estimators L„i (trimmed) and L,„ 2 -, sample sizes 


n= 100 




^n2 

iterration 

mean median minimum 

mean median minimum 

2 

1.001813 1.000857 0.896048 

1.000812 1.000303 0.857210 

3 

1.001827 1.001157 0.824105 

1.000556 1.000109 0.845643 

4 

1.001377 1.000619 0.810082 

1.000360 0.999912 0.844578 

5 

1.000887 0.999840 0.796776 

1.000260 0.999766 0.844182 

6 

1.000372 0.999121 0.777458 

1.000214 0.999723 0.843850 

7 

1.000024 0.999086 0.756777 

1.000196 0.999740 0.843933 

8 

0.999881 0.998921 0.756555 

1.000187 0.999714 0.843928 

9 

0.999806 0.998907 0.756655 

1.000184 0.999726 0.843917 

10 

0.999753 0.998921 0.756655 

1.000183 0.999726 0.843897 
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Table 6. t-distribution: The 25%-quantile, 75%-quantile and max of D-efficiency in 
the sample of5 000 replications of estimators L„i (trimmed) andT,,,!, sample sizes 


n = 100 



*"h1 

^n2 

iterration 

25%-quantile 75%-quantile 

max 

25%-quantile 75%-quantile 

2 

0.980004 

1.023251 

1.157852 

0.986053 

1.014785 

1.194658 

3 

0.972157 

1.030707 

1.212343 

0.983882 

1.016369 

1.212103 

4 

0.967982 

1.032541 

1.214165 

0.983518 

1.016634 

1.214597 

5 

0.967034 

1.032940 

1.222945 

0.983380 

1.016629 

1.215499 

6 

0.966116 

1.032850 

1.260563 

0.983366 

1.016565 

1.215862 

7 

0.965784 

1.032895 

1.261640 

0.983391 

1.016572 

1.216013 

8 

0.965561 

1.032746 

1.261675 

0.983406 

1.016537 

1.216332 

9 

0.965436 

1.032643 

1.261675 

0.983380 

1.016522 

1.216431 

10 

0.965404 

1.032680 

1.261675 

0.983406 

1.016522 

1.216579 


Although resembles the NN-estimator, its behavior for f-distribution reveals its 
robustness no less than . For multivariate normal distribution, both lIP and L^n'^ 
seem to be doing well against outliers. Figures |4] and |5] illustrate this feature in a 
visible way. 
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Fig. 1 Normal distribution: Box-plots of the 5 000 estimated values of 9i(= 1) (top), 62 (= 2) 
(middle) and 83 (= — 1 ) (bottom) for the mean and median. 
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Fig. 2 Normal distribution: Box-plots of the 5 000 estimated values of D-efficiency for the L 
(5) ,(10)^j^(2) ,(5) ^(10) 
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Fig. 3 Normal distribution: Mahalanobis distances represented by co-axial ellipses centered at the 
mean X (black), at the trimmed L„i-estimator (blue) and at the L„ 2 -estimator (red). All simulated 
bivariate data with the every tenth contour are illustrated on the left, detail of the center with the 
first ten contours is on the right. 
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Fig. 4 f-distribution: Box-plots of the 5 000 estimated values of 0i(= 1) (top), 02 (= 2) (middle) 
and 03 (= —1) (bottom) for the L^ 2 ^, L^ 2 ^, Lj^ 2 ^\ mean and median. 
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Fig. 5 ^-distribution: Box-plots of the 5 000 estimated values of D-efficiency for the L 
(10) j(2) ,(5) ^(10) 
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Fig. 6 f-distribution: Mahalanobis distances represented by co-axial ellipses centered at the mean 
X (black), at the trimmed L„i-estimator (blue) and at the L„ 2 -estimator (red). All simulated bivari¬ 
ate data with the every tenth contour are illustrated on the left, detail of the center with the first ten 
contours is on the right. 






















